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Abstract 

We study the large limit of the matrix valued Gross-Neveu model in 
2 < d < 4 dimensions. The method employed is a combination of the ap- 
proximate recursion formula of Polyakov and Wilson with the solution to the 
zero dimensional large counting problem of Makeenko and Zarembo. The 
model is found to have a phase transition at a finite value for the critical 
temperature and the critical exponents are approximated hy v = \/{2{d — 2)) 
and r] = d — 2. We test the validity of the approximation by applying it to the 
usual vector models where it is found to yield exact results to leading order 
in l/N. 
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1 Introduction 



In spite of more than twenty years of efforts, we are still lacking a quantitative 
understanding of the large N limit of matrix valued field theories in more than two 
dimensions. A solution to this problem would be extremely welcome because of its 
twofold application to gauge theory on one side (as pioneered by 't Hooft in and 
string theory and random surfaces on the other (see for a review of the solution 
to lower dimensional systems and for a list of references). 

This state of affairs is to be compared with our successful understanding of vector 
valued models 0. The difference between these two problems is more technical 
than conceptual; for vector models, where the number of fields grows like N, we can 
use saddle point techniques, whereas this is not directly possible if the number of 
fields grows like A^^. This difficulty has been circumvented in dimension d < 1 @] 
and, for gauge theory, = 2 but the naive extension to higher dimensions faces 
formidable computational problems, although recently some progress has been made 
in this direction by pursuing the master field approach [0. 

In a previous paper ^j, we have proposed an approximation to the d dimensional 
problem that combines the known solution to the corresponding zero dimensional 



large N problem and the Polyakov- Wilson recursion formula 0, [jTO[- We have 
applied this approximation to the three dimensional hermitian matrix model and 
found the values for the critical exponents to be 77 = 0.20 and u = 0.67. 

In this paper we extend our investigation to the study of the matrix valued 
Gross-Neveu model in 2 < < 4 dimensions, described by the action 

S = J d'^x tr (^^ + ^^^^"^^ . (1) 



The action has the same form as the one for the original model []TT| with the crucial 
difference that the fields are now arranged into a matrix rather than a column 
vector. This means that the field content is, ignoring the spin indices, ^l{x), \E'^(x), 
k,l = 1, ■ ■ ■ , N). It might be disturbing to some people to have fermionic fields 
living in fractional dimensions but it is possible to make sense of them with the 
usual prescriptions of first performing the calculations for the Dirac algebra and 
then analytically continue in d. Each element of the matrix can be assumed to be 
a two-component complex spinor but their Dirac algebra will never be written out 
explicitly since, as we shall see, the only relevant properties are the anticommuting 
nature of the fields and the existence of a discrete Z2 symmetry preventing the 
occurrence of a bare mass term. This symmetry can be viewed either as "chirality" 
near two dimensions or "parity" near three dimensions: 

d = 2 : x^x ^ ^ 7ch^ ^ -^7ch 

d = 3 : x-^-a;^^^^^ (2) 

Our main results are the following. The fermionic matrix model has the same 
qualitative phase structure as the ordinary vector model, with a phase transition 
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at some critical value for the inverse temperature u* for 2 < < 4, with u* — > 
as 2. The universal quantities at criticality are however quite different, 

showing that the two models are not in the same universality class. In particular, 
to leading order in 1/iV, the "basic" critical exponents for the matrix model within 
our approximation are 

"=2(d^ ^ = ^-2, (3) 

(to be compared with v = l/{d — 2) and r/ = for the vector model). 

As we will argue towards the end of the paper, it might even be that these 
exponents are closer to the actual value than one would generally expect from an 
approximation of this kind. This is related to the existence of a small parameter (the 
thickness of the integration shell in momentum space), under which the recursion 
formula for the fermionic model is better behaved than its bosonic analogue. 

The paper is organized as follows. In section Q we review the nature of the 
approximation and show how it can be applied to various models. We shall try 
to avoid too many repetitions; for more details see 0. In section |^ we apply our 
approximation to the solvable vector model and compare the results with the exact 
solution. We show that in a particular limiting case it is possible to obtain exact 
results for the exponents. A similar phenomenon has been found for the bosonic case 
in []T^ and we shall compare the two cases. Section ^ is devoted to the study of the 
matrix model. The recursion formula is obtained and the critical properties of the 
system are investigated. In particular, the various critical exponents are computed 
for 2 < (i < 4 within this approximation. In section ^ we conclude with a discussion 
of our results and some future projects. The relevant results for the corresponding 
zero dimensional vector and matrix theories are summarized in appendix ^ and ^ 
respectively. 



2 The nature of the approximation. 

The approximation we proposed in p is a combination of the Polyakov- Wilson 
approximate recursion formula [|10| and the solution to the appropriate zero 



dimensional large N counting problem (in this particular case ||T3[)- To perform 
wilsonian renormalization on a (euclidean) field theory one introduces an ultraviolet 
cut-off A, that we shall always set to one by a choice of scale, and then integrates over 
the Fourier components of the fields having momenta p < p < 1 for some p e]0, 1[. 
By rescaling the fields and the momenta so that the new action has the same form 
as the original one, one obtains a recursion relation for the coupling constants of the 
theory describing the critical behavior of the model. 

In many cases the integration of the fast Fourier modes cannot be performed 
exactly. The approximate recursion formula is perhaps the simplest, albeit uncon- 
trolledQ, way to proceed beyond perturbation theory. It amounts to calculating all 

^We shall see later that perhaps, in some special cases, there is a controlling small parameter. 
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Figure 1: The setting sun is one of those diagrams that are not properly counted 
in the original recursion formula because it has an odd number of external lines 
connecting at some vertex. This distinction plays no role in our formulation. 



Feynman graphs by setting the incoming momenta to zero and approximating the 
propagators and the loop integrals by two constants, b and c respectively. It is clear 
that this reduces the integration problem to a zero dimensional counting problem, 
where the dimensionality of space d appears only parametrically through the rescal- 
ing of fields and momenta at the last step of the renormalizationQ. An important 
difference between the fermionic and the bosonic case is that in the former we have to 
allow for b to be complex; the actual phase of b having to be determined a posteriori 
by the requirement that the critical temperature be real positive. 



In its original formulation ||^ [jTO[, this approximation was performed for a single 
component bosonic field One wrote an integral recursion relation for a function 

v{(f)) of a real variable as the approximation to the full effective potential V[^]. 
This formulation has the advantage of keeping track of the renormalization of all 
the ultralocal terms of the kind / d'^x but it has the drawback of ruling out 

wave function renormalization from the onset (see however [0]) and excluding a 
large class of important diagrams such as the "setting sun" depicted in fig. ^ 

The power of combining the recursion formula with the large N limit is that the 
integral for v{(j)) can be performed explicitly (0 is now a vector or matrix valued 
object) reducing the integral recursion relation to an algebraic one that is analytic 
at the gaussian point because of the nice convergence properties of large N integrals. 
All the leading diagrams in are at least counted properly, including the "setting 
sun" and its friends. Wave function renormalization is easy to introduce in this 
scheme; one simply isolates the contributing diagrams, i.e., those diagrams with two 
external legs not joining at the same vertex, and approximates each of them as: 

D(g2) = D{0) + q^D\0) + ■■■ ^ D{0) - ^^^^(O) + • • ■ for bosons 

= D(0)+ ^D'(0) + ---^D(0)-i ^6D(0) + ---forfermions, (4) 

where b is the constant replacing the propagator as before. Replacing the deriva- 
tive by a constant is in the same spirit as the rest of the approximation since the 
derivative "adds" a propagator to the integral of the Feynman diagram. It is also 
the correct substitution on dimensional grounds. 

■^One could also assume that the constant c depends on the dimensionality of space as in 
but it turns out that the universal quantities are independent of this choice, so we shall ignore this 
dependence altogether. 
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In [^] these approximations were applied to the hermitian matrix model 



= \j d'^x tr {{V^f + ro$') + ^ J d'^x tr $^ (5) 

yielding the approximate recursion relations^, (/ = 0, 1, ■ ■ ■) 

ri + 2{l + ri)fr,{g) 



,4-d„. l + /«o(^) 



in terms of three known functions froid): fuoi.9) and fq^^g) (mass, coupling constant 
and wave function renormalization respectively) depending only on the coupling 
constant of the zero dimensional matrix model 0] 

^(0) = tr (^0^ + ^0^) . (7) 

The relation between g and the physical quantities ri and ui was found by analyz- 
ing the topology of the diagrams, (i.e., by counting the relative number of loops, 
propagators and vertices), and it was found that 

/ \ cui 
g{Uhri) 



By studying @ it was possible to estimate the critical exponents for the three 
dimensional hermitian matrix model to be z/ = 0.67 and t] = 0.20. 

There are two main differences between the bosonic case of [§ and the fermionic 
one treated here, irrespective of whether they are vector or matrix models. The first 
difference is that the symmetry (§) prevents a bare mass for the fermionsQ The 
second is that the four fermions interaction is irrelevant at the gaussian point. The 
outcome of this is that the analogue of the pair of equations (|^) is a single recursion 
relation for ui and the approximate renormalization group flow takes place on a line 
rather than on a plane (see fig. ^ for which the gaussian point is attractive. For 
both vector and matrix models, we will show that there is another fixed point u* 
that is repulsive, hence corresponding to a phase transition, and compute the critical 
exponents at that point. In a language perhaps more familiar to particle physicists 
we could call the gaussian point infrared stable and u* ultraviolet stable. There is 
no possibility of confusion as long as we are dealing with a one dimensional flow. 

In [0] a similar calculation was performed for the bosonic vector and matrix 



models where the thickness 1 — p of the integration shell was kept arbitrary. It was 



^In II we had set 6 = 1/(1 + r) and p = 1/2. 

^This symmetry is of course spontaneously broken and the fermions acquire a mass, but all we 
need here is that no Feynman graph gives rise to a term « after integration. 
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found that the critical exponents were mildly dependent on the unphysical value 
of p; an indication that some non universal quantities had "sneaked in" during the 
approximation. By comparison with the bosonic vector model, it was then argued 
that the exponents might become exact in the limit p — >■ 0. We shall see that a 
similar phenomenon occurs here as well, but for the opposite limit p ^ 1. We try 
to give an explanation of this curious fact in section ^ and ||. 



3 Application to the vector model. 

To test our technique and to gain some more insights, let us begin by considering 
the original model [11] with vector valued fermions ^i, \E'"', (i, j = 1, ■ ■ ■ , A^): 

s = jd'x (9) 

This model is very well understood and its critical exponents have been extensively 



studied [15|. To leading order they are rj = and u = l/(d — 2). The discrete chiral 
symmetry is always broken for d = 2, where the model is asymptotically free and 
the fermions are massive. For 2 < d < A, the model has a phase transition at some 
finite value of the temperature, above which chiral symmetry is restored. 

Let us split the field \E' into the sum0 of = + \E'/ carrying only slow (0 < 
p < p) and fast {p < p < 1) Fourier components and define, as usual, 

S[^s + ^f] = S[^s] + ^/] + S[^f]. (10) 

The term (t[^s, ^/] is the part of S that does not factorize under the decomposition, 
namely, 

^[^s, ^f] = ^J d'x {^s^s){^f'^f) + ■ ■ • , (11) 

where the dots represent terms like, for instance, (^'s\E'/)(^s\E'/) that would give 
sub leading contribution in 1/A^ after integrating out the fast modes. This fact 
should be very familiar to the large N experts given the index structure of the 
"cactus" diagrams (see fig. 0). 

We then define a new action S' obtained by integrating out the fast mode^: 

and perform the integral to leading order in 1/A^, keeping the terms 0(\E'^) to all 
orders in perturbation theory, to obtain 

^'[vi>,] = S[^s] + ^Jd''xd'y{^,^,){x){^,^,)iy)x 

( < l^i^f^f){x){^f^f)iy) >conn. - < ^^fix)^fiy) >Ln. ),(13) 

^To keep the notation simple, we indicate only 4" and not ^ but it should always be understood 
that they are independent and subject to the same treatment. 

^The use of the measure exp(S') instead of exp(— 5) is allowed for Grassmann integrals. 
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Figure 2: Among these three simple diagrams, only the first one is leading. The 
thin external lines represent the slow fields, the thick internal ones the fast fields. 
The only vertex allowed between the two in a leading diagram is {'^s'^s){^ f"^ f)- Of 
course, there can be any number of vertices with four internal lines. 



where the connected Green functions are defined with the full measure exp(S'[\E'/]) to 
all orders in uq. Eq. ( p!3[ ) is the exact correction to the four fermions term to leading 
order in 1/N. We claim that if we approximate (0) by making the assumption 
listed in section ^ we obtain: 



(fx 



Mo ^4 



2N 



9o 



(14) 



where r4 is the amputated one particle irreducible four point function for the anal- 
ogous zero dimensional vector model s{ijj) = ipip + g / {2N){%l)'il)Y (see appendix 
and the physical coupling constant uq is related to the zero dimensional one by 



go = uocb'^. 



(15) 



A very explicit way to go from (|T^) to ([T^ is to make the approximation directly 
in (ll^) and then use the zero dimensional Schwinger- Dyson equation 



r4{g)=g + g\CM-Ci{g)). 



(16) 



One can check directly the validity of (|T^ expression by expanding 5" to the 
first few orders in perturbation theory It should be clear however that, after setting 
the external momenta to zero, the only surviving diagrams are the one particle 
irreducible onesQand that the relative powers L, P and V of loops, propagators and 
vertices are controlled by the following topological formulas valid for the four point 
function: 

2V = 2 + P and V - P + L = 1. (17) 

Solving for P and L we see that a diagram with V vertices gives a contribution 
proportional to 



V ^ 0V-1) ^ ^V-l 



v-\ 



as properly encoded in the power series of (0) 

^AU internal propagators vanish for q < p. 
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The final step is to rescale coordinates and fields so that action ([T^) has the same 
cut-off in dimensionless units. Since no wave function renormalization appears to 
leading order for vector models, the rescaling is the "classical" one: x —* p~^x and 
p('^-i)/2^ to get the new coupling constant 

u, = p'-\o X (19) 
90 

It is actually more convenient to multiply both sides by and express everything 
in terms of gi = Uicb'^, {I = 0, 1,2, ■ ■ ■). Thus, the expression for the approximate 
recursion formula for the vector valued Gross-Neveu model becomes 

gi^, = p^-^r,igi). (20) 

The behavior of ( pOD is shown in fig. ^ There are two fixed points; the gaussian 
one and gi = gi+i = g*. What is important is that g* is unstable under gi — > gi+i thus 
defining a phase transition. This is consistent with the fact that the four fermions 
term is irrelevant at the origin and makes the renormalization group flow look as 
shown in fig. ^ for the physical coupling constant u. 

The fact that the critical value of g is negative forces us to choose b to be purely 
imaginary in order to have a positive critical inverse temperature u*. There is no 
danger in doing this as b does not have any direct physical meaning. One can 
solve (|20|) explicitly for the critical point g* and compute the exponent u at that 
point by linearization 0: 

logp logp ^2^) 



A-2d\ 



\og{pd-^T'^(g*)) log(2 - 3p2-d + 2p' 

The exponent u is always positive and independent on the unphysical quantities b 
and c but still depends on the (also unphysical) parameter p. We must admit that 
for, say, p = 1/2, formula ( pTj ) does not do a very good job, giving, for instance, 
u = 1/2 at (i = 3. However, as we let the integration shell become thinner, i.e., 
p 1, formula ( pi] ) becomes exact, i.e., u — > l/{d — 2). This is the analogue 



phenomenon found in ||I2[ for the bosonic case. Its occurrence for vector models is 
not too surprising. For vector models, the leading value for u is obtained by a one 
loop calculation, the "gap equation" . If there is a limiting value for p for which the 
critical coupling u* becomes small, the two calculations are bound to coincide. If 
anything, this observation provides a useful check of our calculation. Also, the fact 



that for bosons [12[ the exact value is reached as p ^ rather that p — 1 is a direct 
expression of the fact that, there, ($^)^ is relevant. 

In the fermionic case, we can think of p as a sort of "control" parameter^ for 
the (generically uncontrolled) approximate recursion formula, in the sense that by 
choosing a very thin integration shell in momentum space, we can keep the critical 



^To be precise, the control parameter is the thickness of the integration shell 1 
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Figure 3: Plot of eq. (^). For clarity, we plot the difference p'^ '^^4(9) ^ 9 ^s. g 
for d = 3 and p = 0.9. The function is real for all g > —1/4; the range shown is 
-0.13 < g < +0.1. 




Figure 4: The qualitative behavior of the renormalization group flow for both vector 
and matrix valued fermionic models. The arrows represent the change in the physical 
coupling ui under iteration / — > / + 1. 
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point near the origin and get exact results for the exponents. Of course, for a thin 
shell the critical point is kept near the origin in the bosonic case as well, but, there, 
the exact exponents are recovered in the opposite limit because of the different 
scaling dimension of the interaction; the interpretation of p as a control parameter 
becomes then more problematic. We now move to the case of matrix valued fields, 
where we shall see that the difference between the bosonic and fermionic case is even 
more significant. 



4 The matrix valued Gross-Neveu model. 

The model of real interest in this paper is described by the action (0) presented in 
the introduction. Since the solution to the corresponding zero dimensional problem 
is known |jl3|], there are no obstacles to extending our technique to this system. 



We can still decompose and \l/ into fast and slow modes, as before, with the 
assumption that all fields are matrix valued. We can still write equation (p!OD, but 
now 

+^f^s^f'^s + 2f /^A*, + 2^',^/'l',^,y (22) 
Although we will never need this in practice, it should be mentioned that some 



of the terms in (22) do not contribute to leading order: the fifth and sixth terms give 
rise to diagrams that are subleading in and the last two terms do not contribute 
to the four point function because they only appear in one particle reducible diagrams 



that vanish by momentum conservation. Hence, only the first four terms in ( p2D 
survive. One could perform the integral over f and find the form of 5" in terms of 
the exact connected Green functions. We skip this tedious calculation because it is 
a repetition of what was done in section | and because a detailed calculation for the 
bosonic case is presented in [§. Instead, we shall argue directly from the diagrams 
what the answer should be. 

We need to compute coupling constant and wave function renormalizationj^. 
Coupling constant renormalization is essentially the same story as for the vector 
model; by setting the external momenta to zero, we select the one particle irreducible 
diagrams, i.e., uq — *■ uoT4{go)/ go where now r4 is the vertex function of the fermionic 
matrix model, computed in appendix The relation go = uqcIP' is unchanged 
because the topological formulas ([TtD are model independent^. 

As for wave function renormalization, we need to isolate the contributing dia- 
grams and apply the approximation (^). We let g and analyze those diagrams 

^"Contrary to the vector model, wave function renormalization appears to leading order in the 
matrix model. 

^^In dl^), the integer L refers to the number of momentum loops, not color loops. 
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Figure 5: The full contribution to wave function renormalization conies from dia- 
gram B. Diagrams like A should be ruled out by symmetry considerations but they 
already vanish in the zero dimensional matrix theory. In this figure, dark circles 
represent the connected Green functions and the light circle the irreducible four 
point vertex. 

that retain some dependence on q in this hmit. Those are all the one particle irre- 
ducible diagrams with two external legs not connected at the same vertex, as shown 
in fig. I^B. All the diagrams that are not one particle irreducible vanish by momen- 
tum conservation because the internal propagators are non zero only for p < g < 1, 
and those diagrams where the two external lines connect at the same vertex (fig. |^) 
do not have any dependence on q. Recalling the symmetry we expect these dia- 
grams to vanish anyway in the exact solution, together with the q independent part 
of all other two point functions. 

When using the zero dimensional results, the symmetry (0) must be imposed 
"by hand"0. However, for the case of the matrix model we encounter a surprise: 
diagram ^ vanishes already in zero dimensions because of the color algebra. This 
fact can be easily checked to one loop in perturbation theory: 



and it is related to another curious property of the fermionic matrix model described 



There is a zero dimensional Z4 discrete symmetry 4' ^ — > — but this is not what we 
want because, for instance, the term tr (^'^'4'^') is not invariant under this symmetry. 



< ^l^i tr ^^^^ > lo = N6l6i{l + 1 - 1 - 1) = 0, 



(23) 
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yielding the integral over 



(25) 



The relative powers of Uq, b and c can be checked as before; the topological 
formulas that apply to the two point function are 

2V=1 + P and V - P + L = 1, (26) 

yielding a contribution for a diagram with V vertices proportional to 
Multiplying by the extra propagator b, to approximate the derivative as in (^), gives 
a contribution to wave function renormalization proportional to u^c^b"^^ = g^. 

The renormalization group recursion relations are obtained as before by rescaling 
X p^^x and p'^'^^^^/^(2 — T2{go))~^^^'^ and redefining the new coupling 

m = p'-'uoX (27) 
^0 (2 -12(^0)) 

Better yet, we multiply ( pTj) by cb"^ to obtain, for I = 0,1 ■ ■ ■, 

9l^, = p'~'-^^^^p'-'R{9l), (28) 
2 - T2{gi) 

where the last equation defines R{g). 

The form of R{g) can be found after a fair amount of tedious algebra to be 

n(n) = 9PiP + ^){P + Q? .... 

"^^y^ + 4)3(2/52 + + 6)' ^ ^ 

where P{g) is the solution to g'^ [3^ — [3"^ — ^(3 — 12 = 0, /3(0) = —2 necessary to satisfy 
the one cut ansatz to the Riemann-Hilbert problem as described in appendix 0. 

By setting gi = gi+i = g* we see that (|^) admits two purely imaginary solutions 
g* = ±ig*. This is fairly easy to show: the function R{g) is odd {f5 is even) and 
real near the origin, where it has a power series R{g) = g — g^ + 0{g^). Hence, for 
p — ^ 1~ we obtain g*"^ ^ 1 — p'^~'^ < leading to two imaginary solutions. This can 
be seen more generally by setting R{ig) = iR{g)- By virtue of R{g) being real and 
odd, R{g) is also real and odd and eq. ( P^ ) becomes 

~gi+i = p^-'m). (30) 

The plot of (^) is shown in fig. ^for the region in which R is real. One can see that 
there are only two Z2 symmetric non trivial solution, both of which are unstable. 
The requirement that u* be real positive forces us to fix b^ to be purely imaginary 
and to discard one of the two solution as unphysical. The renormalization group 
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flow involving the gaussian point and u* is qualitatively the same as the one depicted 
in fig. |. 

There are other solutions to (|2^) in the whole complex plane but the above 
solution is the only one that can be really trusted by virtue of reducing to the 
gaussian when the integration shell becomes thinner (p — >■ 1). In fact we see that, 
just as for the fermionic vector model in section ^, the parameter p provides a way 
of controlling the approximation by keeping the critical point near the gaussian 
one for thin integration shells. The situation is more favorable than in the bosonic 
case |]12[ because it is in this region that the exponents for the fermionic vector model 
become exact. It is possible to linearize ( pH]) numerically and extract the exponent 
V. However, as we have argued that the results should be reliable for p 1, it is 
enough to consider the first non trivial term in the Taylor expansion near the origin: 

gi+.^p^-^gi-g^). (31) 

To illustrate why we expect a universal behavior different from the one of the vector 
model, let us recall that in the vector case the Tailor expansion of (^) contains a 
quadratic term: (yf/+i ^ p'^~'^{gi — gf)- In the most general case 

gi+,^p''-\gi + Ag'^), (32) 

one obtains 

^ ~ — ; -, , , ^, -, z-r-, r as p — i> 1. (33) 

log(n- (n- l)p'^-2) (n-l)(rf-2) ^ ^ ' 

For n = 2 we recover the exponent for the vector model and for n = 3 we prove our 
first claim on the matrix model: v = l/{2{d — 2)). 

In a similar way, it is possible to evaluate the wave function renormalization at 
the critical point and to extract the exponent 

iog(2 - r2(^*)) , _ ^ 

T] = > a — 2 as p — > 1. (34) 

logp 

If we are willing to be optimistic, we can use the combined values of u and r] to 
compute the other magnetic exponents from the well known scaling relations: 

3c? -8 1 A-d 2 

a = —— -, p = -, 7 = — and o = . (35) 

2(d-2)' 2' ' 2{d-2) d-2 ^ ^ 

It should be possible to test the validity of some of these results by a computer 
simulation. 



5 Conclusions. 

We have shown how to compute the critical exponents for the matrix valued Gross- 
Neveu model by using a combination of the approximate recursion formula and 
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the zero dimensional large N limit. The success of this technique in predicting the 
exponents for the ordinary vector models (both bosonic and fermionic) gives us more 
confidence in the approximation. The recursion relation for the fermionic models is 
even better behaved that for the bosonic ones; the thickness of the integration shell 
controls the distance of the critical point from zero. 

Since we have used a Taylor series expansion to obtain the values of the criti- 
cal exponents, one might be tempted to say that a perturbative calculation should 
suffice. This is not the case because the existence of a non trivial fixed point can 
never be established by perturbation theory alone; one must rely on some non per- 
turbative treatment such as the one used in this paper. However, it is true that for 
a certain range of the parameter p, the fixed point falls near the origin and then 
there is no harm in doing a Taylor expansion. In a sense, the situation is similar 
that in the e expansion |jl6| where the controlling parameter e can be chosen in such 
a way that the fixed point also falls near the origin, thus allowing one to expand 
in powers of the coupling constant. The parameter p plays a role similar to e but 
only for fermionic models. One crucial difference is that e is a physical parameter 
(essentially, the dimension of space) whereas p is not. It would be interesting to 
develop a self consistent way to show that is is possible to neglect higher powers of 
from the exact renormalization group equations. 

Two directions worth exploring in the future are the connection with randomly 
triangulated surfaces and the application of similar methods to gauge theories. In 
both cases, the simple approximations used in this paper need to be modified for 
many reasons; the critical points studied here and in ^ are not the ones correspond- 
ing to the continuum limit of the string partition function and gauge symmetry is 
not preserved in the presence of a sharp cut-off. A solution to either of these prob- 
lems, allowing one to carry out the spirit of the approximation for these systems, 
would be extremely welcome. 



Acknowledgements. 

We would like to thank J. Ambj0rn, P. Damgaard, S. Nishigaki and G. Semenoff for 
discussions. 



A The zero dimensional fermionic vector model. 

In this appendix we quickly review the relevant results for the fermionic vector 
model. We will be brief but try to be self contained; for more details see |T^. The 
great simplification for the vector models is that their Schwinger-Dyson equations 
can be solved exactly in a closed form, contrary to the matrix model case where, 
even after imposing the Schwinger-Dyson equations, one is still left with some un- 
determined correlators. Hence, for the vector model, we can obtain the quantities 
we need without doing any real work! 
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For any quantity O we define 

^ ^~ J dip dip e^^+^^^^^^ 
It is a simple matter to check that the first two Schwinger-Dyson equations read 

< j^ipip >+g< ■^{^i'f > -1 = 

< ^(^^)' > +^7 < ^(V'^f > - (l - ^) < > = 0. (37) 

Eq. ( p?!) is exact, i.e., vahd for all finite A^'s. Now we define the connected correlation 
functions by cluster decomposition: 



< ip'ipj > = < ip'ipj > 



J ^conn. 



< ip^ipjip^ipi > = < ip'ipj >conn.< ip'^'ipi >conn. " cross term 

+ < Ip'lpjlp^lpl > 

conn. 

< Ip'lpjlp^lpilp'^lp^ > = < Ip'ipj >conn.< >conn.< ^'"'V^n >conn. ± flVC C.t. 

+ < Ip'ipj >conn.< Ip'^lpllp'^lpn >conn. ± eight C.t. 

+ < iP'lPjlP^lPilP'^lPn >conn. • (38) 

For our purposes, we need to keep track of terms of order in (|38|) and this can 
be done by introducing three quantities 6*2(5'), 6*2 (5) and C^i^g) through 

<#^, >conn. = 5]{c,{g) + ^C,{g)^+0{j^) 

< ^V#V^^ >conn. = {5]5l-5\5')l^C,{g)+0{j^) 
<lP'lPjlP%lP"'lPn>conn. = O(^). (39) 

Summing over the color indices in (|38| ) and substituting into (|37|) we obtain, after 
keeping careful track of all powers of 

C2 + gCl- 1 = 

C2 + g{2C2C2 + C4 - C|) = 

2C2C2 + C4 - Cl + ^gC2 {C2C2 + Ci-Cl)-C2 + C2 = 0. (40) 

Eq. ( ^0|) forms a closed system of equations that can be solved for the quantities 
that we need: 

CM = ^^^^ 

2g 



{l + 2g-./T+Tgf 
^'^'^ = 2,^(1 + 4, -VI + 4,) - ^''^ 
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By a decomposition similar to (^) it can be immediately checked that the relations 
between the connected Green functions C2 and C4 and the vertex functions r2 and 
r4 are the usual ones: 

r,(,) ^ ^ i±4±3« 

Ug) = cMCM-' = i(i + -^y (42) 

where the only unusual thing is, perhaps, the relative sign between C4 and r4 that 
is dictated by the choice of the measure exp(+S'). The vertex function r4 is the one 
used in the computation for the vector model in section ^ 



B The zero dimensional fermionic matrix model. 

We now move on to the fermionic matrix model for which we have to do a bit more 
work since the Schwinger-Dyson equations do not close. For any quantity O we now 
define 

< C >= ^ ^ ^ _ , - - . . 43 

The only Schwinger-Dyson equation we shall need is 

< ^ tr (^^) >+g<^tT (xp^tptp) > -1 = 0. (44) 

The trick of cluster decomposition of appendix |^ will not work here because the 
connected components are of the same order as the disconnected ones. Instead, we 
follow and define a density^ p(/i) 



1 /" 
< tr >= J dfi p{fi)ij\ (45) 

and a generating function u{z) 

.(.)^<it.i^]>^[,,I^. (46) 



N J z - a 



We shall not repeat the argument here but it is well known |]T3[ that the density 
p{fi) is the same as the density for the Penner model i.e., it has support on a 
finite complex arch on which it satisfies 



^•^ A matrix with anticommuting elements cannot be diagonalizcd and p no longer has the meaning 
of eigenvalue density. 
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The one cut solution to the Riemann-Hilbert problem (^) is: 

^.) = i-i-f +(f-^)y?T^, (48) 

where a and /3 are determined by the asymptotic condition Ci;(z) — at infinity. 
The quantity /3 is the same as the one used in section ^. After some tedious algebra, 
it is possible to eliminate a from the problem and to write the two point function 
as 

<^-(«>-^' 

where (5 is the solution to (7^/3'* — — 8/5 — 12 = with the initial condition 
/3 = — 2 at (7 = necessary to ensure the proper gaussian limit. The four point 
function can be read off from (0). 

To conclude, we wish to show how all this is related to the connected Green 
functions and C\ and the vertex functions r2 and This is achieved by simple 
cluster decomposition: 

< ^j^f > = < ^)^\ >conn. 
< V^iV^fCV'? > = < ^i^f >conn.< C< >conn. 

- < >conn.< >conn. 

+ < ^;-^f C< >conn. (50) 



and by the usual definition 



1^ 

N' 



Substituting (0) into (^) and taking the trace, we see that the disconnected piece 
drops out because of the antisymmetry in the color indices yielding 

CM ^ <-tr(^^)>=— ^ 

1 - - /?2 -I- -I- fi 

^4(^?) - <^tr(^^^^)>= ^^+^+ . (52) 

The result for C4 is rather curious and applies only to fermionic matrix models; for 
all other models a disconnected term C\ would survive after tracing. 

Finally, by a decomposition similar to (0), we see that the vertex functions are 
again related to and C4 as in appendix 0: 

VM = CM-^ and YM = CMCM-\ (53) 

These are the quantities used in section ^ to obtain the renormalization group equa- 
tions. 
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